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Abstract 

In the framework of the fcT^-factorization approach, the production of unpolarized Drell- 
Yan lepton pair at high energies is studied. The consideration is based on the 0{a) and 
0{aas) off-shell partonic matrix elements with virtual photon 7* and Z boson exchange. 
The calculations include leptonic decays of Z bosons with full spin correlations as well as 
7* — Z interference. The unintegrated parton densities in a proton are determined by the 
Kimber-Martin-Ryskin prescription. Our numerical predictions are compared with the data 
taken by the D0, CDF and CMS collaborations at the Tevatron and LHC energies. Special 
attention is put on the specific angular distributions measured very recently by the CDF 
collaboration for the first time. 

PACS number(s): 12.38.-t, 12.15.Ji 



1 Introduction 

With the start of the LHC experiments, high energy particle physics entered a new era. 
The LHC opens a new kinematic regime where the number of novel physical phenomena 
can occur. One of the most useful tools to study hadronic interactions at high energies is 
the Drell-Yan dilepton production, in which quark-antiquark annihilation form intermediate 
virtual photon or Z bosons decaying to lepton pairs. This process is presently of considerable 
interest from both experimental and theoretical points of view. In particular, Drell-Yan pair 
production is a unique process which offers high sensitivity to the parton (quark and gluon) 
distributions in a proton. It provides a major source of background to a number of processes 
such as Higgs, tt pair, di-boson or W and Z' bosons production (and other processes beyond 
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the SM) studied at hadron colliders. Dilepton production has a large cross section and clean 
signature in the detectors and therefore it is used for monitoring of the collider luminosity and 
calibration of detectors. Moreover, it is an important reference process for measurements of 
electroweak boson properties at hadron colliders. Therefore it is essential to have an accurate 
QCD predictions for corresponding cross sections and related kinematical distributions. 

Theoretical investigations of Drell-Yan pair production have own long story. It is one of 
the few processes in hadron-hadron collisions where the collinear QCD factorization has been 
rigorously proven [IHl]- Within this framework, the NLO pQCD calculations of inclusive 
cross sections have been performed [5H7], and later it was done on up to NNLO accuracy [8|l9] . 
Recently fully exclusive NNLO pQCD calculations became available, including the leptonic 
decay of intermediate Z boson [TOHTS] . The results of these calculations agree with the 
Tevatron and LHC data within the theoretical and experimental uncertainties. Of course, 
typically for collinear QCD factorization, perturbative calculations diverge at small dilepton 
transverse momenta pr ^ ^ (where M is the invariant mass of produced lepton pair) with 
terms proportional to InM/p^ appearing due to soft and collinear gluon emission. Therefore 
special soft gluon resummation technique p^] - |20] should be used to make QCD predictions at 
low pt- Such soft gluon resummation can be performed either in the transverse momentum 
space [21] or in the Fourier conjugate impact parameter space [22] • Differences between 
the two formalisms are discussed in [23]. The traditional calculations combine fixed-order 
perturbation theory with analytic resummation and some matching criterion. 

An alternative description can be achieved within the framework of the /cT-factorization 
approach of QCD [24J. This approach is based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
[25j or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [26] equations and provides solid the- 
oretical grounds for the effects of initial gluon radiation and intrinsic parton transverse 
momentum fcy. A detailed description and discussion of the /cr-factorization formalism can 
be found, for example, in reviews [2Z|. Here we only mention that, in contrast with the 
collinear approximation of QCD, the initial gluon emissions in the /cT-factorization approach 
generate the finite dilepton transverse momentum already at Born level. Moreover, the 
soft gluon resummation formulas are the result of the approximate treatment of the solutions 
of CCFM equation, as it was shown in |28j . 

In the present note we apply the /cT-factorization approach to unpolarized Drell-Yan pair 
production in and collisions. A non-collinear factorization theorem for this process has 
been proven [29j for M . Below we assume it in a wide range of pr for phenomenological 

purpose^. First application of /cT^-factorization approach to lepton pair production has 
been performed in [31], where authors have considered only diagrams with virtual photon 
exchange and concentrated mostly on the rather low energies covered by the RHIC and UAl 
experiments. More general consideration of high energy resummation for Drell-Yan processes 
was done in [3^. Our main goal is to give a systematic analysis of the Tevatron data [33n38j 
and first LHC measurements [SH] performed by the CMS collaboration. The consideration 
is based on the 0[oi) and 0[aa^ off-shell (depending on the non-zero transverse momenta 
of incoming partons) production amplitudes where we take into account both 7* and Z 
boson exchange. Our calculations include also leptonic decays of Z bosons with full spin 
correlations. Thus we easily can produce various kinematical distribution and apply cuts, 

^For discussion of the fcT-factorization for high energy processes see, for example, [30] . 
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analogous to experimental ones. Specially we study the angular distributions of produced 
lepton pair measured very recently by the CDF collaboration for the first time and 
investigate the different sources of theoretical uncertainties. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the fc]-- factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 



2 Theoretical framework 



There are three subprocesses which describe Drell-Yan pair production at order of 0{a) 
and 0{aas): 

q + q^Y/Z ^l^ + 1- (1) 
q + 9* ^1*/Z + q^l+ + 1- + q (2) 
q + q^Y/Z + g^l+ + l-+g (3) 

The corresponding Feynman graphs are shown on the Fig. 1. Note that in the framework of 
/cT-factorization approach contribution from the subprocess (3) is already taken into account 
by the quark-antiquark annihilation (1) due to the initial state gluon radiation. Therefore 
below to avoid the double counting we consider only the subprocesses (1) and (2). It is in 
a contrast with the coUinear QCD factorization where contributions from the subprocesses 
(1) — (3) should be taken into account separately. 

Let us start from the kinematics. We denote the 4-momenta of incoming partons and 
outgoing leptons by ki, k2, pi and p2- The initial hadrons have the 4-momenta p^^^ and p^'^\ 
and the final state quark in (2) has 4-momentum p^. In the center-of-mass frame of colliding 
particles we can write: 

P^^) = ^(1,0,0,1), p(2) = ^(1,0,0,-1) (4) 

where ^/s is the total energy of the process under consideration and we neglect the masses 
of the incoming protons. The initial parton four-momenta in the high energy limit can be 
written as follows: 

ki = xip^^^ + kiT, k2 = XiP^'^^ + fc2T, (5) 

where kyr and k2T are the corresponding transverse 4-momenta. It is important that k^^" — 
~^iT 7^ 0) = —k\Y 7^ 0. From the conservation laws we can easily obtain the following 
relations for annihilation subprocess (1): 

kiT + = Pit + P2T, (6) 

xia/s = miye^^ m2Te^^ (7) 
X2\rs = miTe~y^ + m2Te~^^ (8) 
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and the similar ones for QCD Compton subprocess (2): 



kiT + kar = Pit + P2T + Psr, (9) 

Xl^/s = mire^^ + mare^^ + ^tistc^^, (10) 
X2\/s = irtiTe'^^ + m2Te"^^ + m3Te"^^ (11) 



where piT, P2T and psr are the transverse momenta of produced particles, ?/2 and ?/3 
are their center-of-mass rapidities and miT, '^2T and msr are the corresponding transverse 
masses, i.e. mfrp = mf + 'plj,. The matrix elements of (1) and (2) can be presented as follows: 

Ml = ie'^eqVs,{k2)YUs2{kl) — UrMYVr2{P2), (12) 

s 

4 COS 



X 



vp?y ) is — vn?y — imzTz) 



(13) 



^2 



(14) 



4 COS"' Ow 

^ . _ {Pl+P2)p{Pl+P2)u \ UrAPl)Y{C^ - C\Y)Vr2ij>2) 

rn% ) iPi + P'iY irrizTz 

where e and Cq are the electron and quark (fractional) electric charges, s = (/ci + A;2)^, dw and 
are the weak and strong charges, mz and Tz are the mass and full decay width of Z boson, 
6\y is the Weinberg mixing angle, and a are the polarization 4-vector and eight-fold color 
index of incoming off-shell gluon, Cy and Ca are the vector and axial constants. Here we 
neglected the masses and virtualities of incoming quarks and took propagator of intermediate 
Z boson in a Breit-Wigner form to avoid an artificial singularities in numerical calculations. 
When we calculate the matrix elements squared, the summation over the incoming off-shell 
gluon polarizations is carried out in according to the /cj^-factorization prescription 



e'^e — k.2rpk.2rp/]<.2rp. (16) 

In the collinear limit, where \\<.2t\ 0, this expression converges to the ordinary ^e'^e*'^ = 
after averaging on the azimuthal angle. In all other respects the evaluation follows the 
standard QCD Feynman rules. The calculation of traces in (12) — (15) is straightforward 
and was done using the algebraic manipulation systems Form |10]. We do not list here 
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the obvious expressions because of lack of space. The obtained expression for Compton 
subprocess (2) coincides with the result [32] . 

To calculate the cross section of Drell-Yan lepton pair production, in according to the /c^- 
factorization theorem, one should convolute off-shell partonic cross sections with the relevant 
unintegrated quark and/or gluon distributions in a proton: 



(17) 



where a*j{xi,X2, k^j., k^j^) is the off-shell partonic cross section and fi(x, k^, /i^) is the unin- 
tegrated parton densities in a proton. The contributions to the total Drell-yan cross section 
from quark- ant iquark annihilation and QCD Compton subprocesses can be easily rewritten 
as follows: 

q 

'\ d(t)2 



a 



:\Ml 



X 



167r(a;ia;2s)^ 

xf\{xi, kiy, i/)fq{x2, kaj., iJ,'^)dplTd\4j.dklj.dyidy2- 



27T 27r 



a 



2567r3(xiX2s) 



■\M 



X 



x/g(xi, kjj^, /x^)/g(x2, kjj,, jj'^)dplj,dplj.dklj.dklj.dyidy2dy: 



d(f)i d(j)2 dipi dip2 



(19) 



2tx 27r 27r 27r 

where 0i, 02, 4^1 and V'2 are the azimuthal angles of initial partons and produced leptons, 
respectively. If we average these expressions over 0i and 02 and take the hmit |kiT| — )■ and 
|k2r| 0, then we recover the corresponding formulas in the coUinear QCD factorization. 

Concerning the unintegrated quark and gluon densities in a proton, we apply the Kimber- 
Martin-Ryskin (KMR) approach |?T] to calculate them. The KMR approach is the formalism 
to construct the unintegrated parton distributions from the known conventional ones. In this 
approximation the unintegrated quark and gluon distributions are given by 



fq{x, k|i, /i^) — Tg(ky, /x^ 



a. 



27r 



X 



X 



dz 



Pqqiz)-q (-, k|) e (A - ^) + P,,iz)-g 
. z \z / z 



(20) 



X 



/g(x, ky,/i^) — Tg(k^,/i^ 



(k^) 



2tt 



X 



X 



± 

I dz 5^P,,(.)|gg,k|)+P,,(.)|^(|,k|)e(A-.) 



(21) 



where Pab{z) are the usual unregulated LO DGLAP splitting functions. The theta functions 
which appears in (20) and (21) imply the angular-ordering constraint A = /i/(/i + jk^^l) 
specifically to the last evolution step to regulate the soft gluon singularities. Numerically, 
for the input we have used leading-order parton densities xq{x, iJ?) and xg{x, iJ?) from recent 
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MSTW'2008 set [42j. The Sudakov form factors Tq{\4, /j,^) and Tg{\4, /j,^) enable us to 
include logarithmic loop corrections to the calculated cross sections. To take into account the 
non-logarithmic loop corrections we use the approach proposed in |33] . It was demonstrated 
that main part of the non-logarithmic loop corrections to the quark-antiquark annihilation 
cross section (1) can be absorbed in the effective i^'-factor: 



K 



exp 



27r 



-IT 



(22) 



where color factor Cp = A/3. A particular choice /x^ = p^^M"^^^ has been proposed 
to eliminate sub-leading logarithmic terms. We choose this scale to evaluate the strong 
coupling constant in (22). 

The multidimensional integrations in (18) and (19) have been performed by the means of 
Monte Carlo technique, using the routine VEGAS [H]. The full C++ code is available from 
the author on request^]. 



3 Numerical results 



We now are in a position to present our numerical results. First we describe our input 
and the kinematic conditions. After we fixed the unintegrated gluon distributions, the cross 
sections (18) and (19) depend on the renormalization and factorization scales fiji and fip. 
Numerically, we set them to be equal to nn = hf = ^M. To estimate the scale uncertainties 
of our calculations we vary the parameter ^ between 1/2 and 2 about the default value ^ = 1. 
Following to I45j, we set mz = 91.1876 GeV, Tz = 2.4952 GeV, sin^ Ow = 0.23122 and use 
the LO formula for the strong coupling constant agifi"^) with = 4 active quark flavors at 
Aqcd = 200 MeV, so that as(M|) = 0.1232. 

The results of our calculations are presented in Figs. 2 — 4 in comparison with the 
D0 [38], CDF [33H37] and CMS data [39]. Solid histograms are obtained by flxing both 
the factorization and renormalization scales at the default value /i = M, whereas the upper 
and lower dashed histograms correspond to the scale variation as it was described above. 
The predicted total cross sections are listed in Table 1. One can see that the Tevatron and 
LHC experimental data are reasonable well described by the /cT-factorization approach in 
the whole range of invariant masses. Our predictions tend to only slightly overestimate the 
rapidity distribution of dilepton pair in the region of Z boson peak 66 < M < 116 GeV, but 
agree with the data within the uncertainties. Specially we point out a good description of 
dilepton transverse momentum distributions measured by the CDF collaboration since this 
observable strongly depends on the unintegrated parton density used. 

The relative contributions of quark-antiquark annihilation and QCD Compton subpro- 
cesses to the Drell-Yan cross sections at the Tevatron and LHC energies are shown in Fig. 5 
as a function of azimuthal angle difference between the transverse momenta of produced lep- 
tons. Note that this observable is singular in the collinear QCD approximation at LO due to 
back-to-back kinematics. It is in a contrast with the fcr-factorization approach, where, as it 

^ lipatov@ theory, sinp . msu.ru 
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Source 


a(66 < M < 116 GeV), pb 


a{M > 116 GeV), pb 


/cT-factorization (KMR) 
NNLO pQCD 
CDF data [36] 


285 ± 31 
227 ±9 
250 ±4 (stat.) ±10 (syst.) 


3.7 ±0.4 
3.3 

4.0 ±0.4 (stat.±syst.) ±0.2 (lumi.) 



Table 1: Total cross sections of Drell-Yan pair production in pp collisions at y/s = 1800 GeV. 



was mentioned above, the finite transverse momentum of dilepton pair is generated already 
in LO quark- ant iquark annihilation (1). We find that latter dominates at high A0 ~ vr for 
both the Tevatron and LHC energies, whereas at A0 < 7r/2 quark-antiquark annihilation 
and QCD Compton subprocesses contribute equally. Note that here we applied no cuts on 
the final-state phase space. 

Now we turn to more detailed analysis of angular distributions in dilepton production. 
The general expression can be described by the polar 6 and azimuthal (p angles of produced 
particles in the dilepton rest frame. When integrated over cos 6' or 0, respectively, the angular 
distribution can be presented as follows: 



~ (l + cos^e) + -Ao(l-3cos^9) + A4Cose, 
a cos 9 2 



da 

— — ~ 1 



/Ss cos 4> + /32 cos : 



(23) 



(24) 



where /^s = SttA^/IQ and f32 = ^2/4. Note that the angular coefficients Aq and A2 are the 
same for 7* or Z boson exchange, and A^ and originate from the 7* — Z interference. 
The Lam- Tung relation [46] Aq = A2 is valid for both quark-antiquark annihilation and 
QCD Compton subprocesses at 0{aas) order. Higher-order QCD calculations [HlllH] as 
well as QCD resummation up to all orders [19] indicate that violations of the Lam- Tung 
relation are small. Very recently the CDF collaboration reported [37] the first measurement 
of the angular coefficients Ao,A2, A^ and A4 in the Z peak region (66 < M < 116 GeV) at 
\/s = 1960 GeV. Below we estimate these coefficients regarding the CDF measurements. Our 
evaluation generally followed the experimental procedure. We have collected the simulated 
events in the specified bins of dilepton transverse momentum, generated the decay lepton 
angular distributions according to the matrix elements (12) — (15) and then applied a two- 
parametric fit based on (23) and (24). The estimated values of angular coefficients in the 
Collins-Soper frame are shown in Fig. 6. We find that our predictions agree well with the 
CDF data as well as collinear QCD predictions listed in [37]. We would like to only remark 
that the latter predict a flat behaviour of ^3 in a whole px range whereas CDF data tends 
to support our predictions (slight decreasing of A^ when we move to large px values). 



7 



Finally, we can conclude that /cr-factorization predictions in general are rather similar 
to ones based on the collinear QCD factorization with the NNLO accuracy. It demonstrates 
again that the /cT^-factorization approach at LO level automatically incorporates a large piece 
of the standard (collinear) high-order corrections |2Z]- It is important for further studies of 
small-x physics at hadron colliders, and, in particular, for searches of effects of new physics 
beyond the SM at the LHC. 

4 Conclusions 

We have investigated unpolarized Drell-Yan lepton pair production in pp and pp collisions 
at the Tevatron and LHC energies within the framework of the /cT-factorization approach. 
Our consideration is based on the 0{a) and 0{aas) off-shell production amplitudes where 
7* and Z boson exchange is taken into account. The calculations include leptonic decays of 
Z bosons with full spin correlations and 7* — Z interference. The unintegrated parton den- 
sities in a proton are determined by the Kimber-Martin-Ryskin prescription. We obtained a 
reasonable well agreement (at a similar level as in the NNLO pQCD) between our predictions 
and the available data taken by the D0, CDF and CMS collaborations. Specially we studied 
the specific angular distributions measured very recently by the CDF collaboration for the 
first time. 
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Figure 2: The total and differential cross sections of the Drell-Yan pair production in pp and 
pp collisions at the Tevatron and LHC as a function of dilepton invariant mass M. The solid 
histograms correspond to the results obtained with the KMR parton densities. The upper 
and lower dashed histograms correspond to scale variations, as it is described in the text. 
The experimental data are from D0 p], CDF [33l|35] and CMS |39]. 




Figure 3: The differential cross sections da/dy of dilepton production at i/i = 1800 TeV 
compared to the CDF data [36]. Notation of all histograms is the same as in Fig. 1. 
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Figure 4: The differential cross sections da/dpT of dilepton production at ^/s = 1800 TeV 
compared to the CDF data [3l]. Notation of all histograms is the same as in Fig. 1. 
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Figure 5: Different contributions to the Drell-Yan pair cross section in pp and pp collisions 
at the Tevatron and LHC energies as a function of produced lepton azimuthal angle differ- 
ence. The dashed and dash-dotted histograms correspond to the contributions from quark- 
antiquark annihilation and QCD Compton subprocesses, respectively. The solid histograms 
represent the sum of these components. 
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Figure 6: Angular coefficients Aq, A2, and of dilepton production as a function of pt 
compared to the CDF data [37]. Solid and two dashed histograms represent fitted values of 
angular coefficients and corresponding uncertainties of fitting procedure. The default scale 
H = M has been applied. 
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